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ABSTRACT: The Heat source/sink and suction/injection effects are studied during the Heat and Mass 
transfer through copper, water nano -fluid along an inclined permeable oscillating flat sheet. The 
governing equations are solved and the influence of various parameters is analyzed. The Rate of heat 
transfer for volume fraction against heat source is also analyzed. 
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List of symbols: 

Bf\ Constant applied magnetic field (Wb m" ) 
C p Specific heat at constant pressure (J kg" 1 K" 1 ) 
Q Gravity acceleration (m s" ) 
/ Current density 

M Dimensionless magnetic field parameter 

71 Dimensionless frequency 

Nu Local Nusselt number 

Pr Prandtl number 

Q Dimensional heat source (kJ s -1 ) 

Dimensionless heat source parameter (kJ s" 1 ) 
Dimensionless suction parameter 
Dimensionless time (s) 
Local temperature of the nano-fluid (K) 
Wall temperature (K) 

Temperature of the ambient nano-fluid (K) 
"Uj w Dimensionless velocity components (m s" 1 ) 
£/(j Characteristic velocity (m s" 1 ) 
w$ Mass flux velocity 

k Thermal conductivity 
Df Diffusivity of water 

D. z Diffusivity of copper 



£ 



Greek symbols: 

& Thermal diffusivity (m s" 1 ) 

Thermal expansion coefficient (K" ) 



£ 



Molecular expansion coefficient 
Dimensionless small quantity (« 1) 

Solid volume fraction of the nano-particles 

Dynamic viscosity (Pa s) 
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2 1 

y/ Kinematic viscosity (m s" ) 
6 Dimensionless temperature 

2 1 

C Electrical conductivity (m s" ) 

a x Stefan-Boltzmann constant 

8 Mean absorption Coefficient 
y Inclination angle of the plate 
p Density 

C Dimensionless diffusion 
Superscript 

- Dimensional quantities 

Subscripts 
f Fluid 
s Solid 
nf Nano-fluid 

I. INTRODUCTION 

Research in the field of Heat and Mass transfer challenging the cooling of many systems used in day to 
day life of mankind.The heat and mass transfer enhances enormously when nano -particles are suspended in 
liquids like water, ethylene glycol etc. This has substantiated by Das, Choi and Patel (2006) in their review 
paper. In this scenario cooling systems demand the very low heat and mass transfer rate through nano - fluids 
and heat mass energy systems like automobiles demanding the high heat and mass transfer rate through nano - 
fluids. 

Kuznetsov and Nield (2010) studied the classical problem of free convection boundary layer flow of a 
viscous and incompressible fluid (Newtonian fluid) past a vertical flat plate to the case of nano -fluids. In these 
papers the authors have used the nano-fluid model proposed by Buongiorno (2006). Although this author 
discovered that seven slip mechanisms take place in the convective transport in nano -fluids, it is only the 
Brownian diffusion and the thermophoresis that are the most important when the turbulent flow effects are 
absent. More recently, Khan and Aziz (2011) studied Natural convection flow of a nano-fluid over a vertical 
plate with uniform surface heat flux. M. A. A. Hamad and I. Pop (2011) presented in their recent paper that the 
solid volume and heat source enhances the heat and mass transfer rate. This brief survey clearly indicates that a 
definitive conclusion regarding the role of nano -particles in enhancing natural convective transport is yet to be 
reached. Recently Anwar et.al (2012) studied the conjugate effects of heat and mass transfer of nano-fluids over 
a non-linear stretching sheet. 

In this paper we aim to investigate the MHD Cu - water nano-fluid flow and the heat and mass transfer 
past a vertical infinite permeable inclined oscillating flat plate under heat and mass source and suction. 

II. MATHEMATICAL FORMULATION 

Consider the unsteady three dimensional free convection flow of a nano-fluid past a vertical permeable 
semi-infinite plate in the presence of an applied magnetic field with constant heat source. We consider a 

Cartesian coordinate system (x,y,z) .The flow is assumed to be in the x direction, which is taken along the 

plate, and z - axis is normal to the plate. We assume that the plate has an oscillatory movement on time t and 

frequency n with the velocity u(0,t), which is given u(0,t) = U 0 (1 + s cos(nt)), where s is a small constant 
parameter ( s « 1) and U 0 is the characteristic velocity. We consider that initially (t < 0) the fluid as well as the 

plate is at rest. A uniform external magnetic field B 0 is taken to be acting along the z -axis. We consider the case 
of a short circuit problem in which the applied electric field E = 0, and also assume that the induced magnetic 
field is small compared to the external magnetic field B 0 . The surface temperature is assumed to have the 
constant value T w while the ambient temperature has the constant value T M , where T w > T^. The conservation 
equation of current density V./=0 gives J z = constant. Since the plate is electrically non-conducting, this 

constant is zero. It is assumed that the plate is infinite in extent and hence all physical quantities do not depend 
on x and y but depend only on z and t , 
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i.e. 



d u d v 

+ — = 0 



dx dy 

It is further assumed that the regular fluid and the suspended nano -particles are in thermal equilibrium 
and no slip occurs between them. Under Bossinesq and boundary layer approximations, the boundary layer 
equations governing the flow and temperature are, 



d w 

iz 

d u 
d T 

d c 
71 



= 0 



•(1) 
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+ w 



+ w 



9 u 
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p nf 



o u 



d z 
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u 



= a 



d 2 T 

0 z 



Q 



•(3) 



= D 



o c 
0 z 



•(4) 



The appropriate initial and boundary conditions for the problem are given by 



u(zj) = 0, T = ,c = for KOVz 



00 



u(0,t) = Uq 



-t . £ , int , — intx 
1h — (e +e ) 

2 



T(0,t) = T w ,c(0,t)=c 



w 



t>0 



w(oo,f)->o,r(oo,r)-^r 00 , c(oo,r) Cao , s « l 

Thermo-Physical properties are related as follows: 



Pnf = Q-t)Pf + tP s > Pnf 



(1-0 



2.5 



, a 



nf 



nf 



(Pc n ) 



p'nf 



k + 2k r — 2S(k r —k ) 



nf 



= k 



f 



k +2kr + (b(kr—k) 
s f r f s J 

The thermo-physical properties (values) of the materials used are as follows 



(6) 



Table I 



•(2) 



.(5) 



Physical Properties 


Water 


Copper(Cu) 


C p (J /kg K) 


4,179 


385 


p (kg In?) 


997.1 


8,933 


k (W/mK) 


0.613 


400 


P T xlO~5 (\/K) 


21 


1.67 


8 x 10 6 (m 2 //i) 
c 


298.2 


3.05 



We consider the solution of Esq. (1) as w = - w 0 (7) 

Where the constant w 0 represents the normal velocity at the plate which is positive for suction (w 0 >0) and 
negative for blowing or injection (w 0 <0).Thus, we introduce the following dimensionless variables: 
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u, 



0 



z=(-*-)Z,t = (-^-)t*, n = (-^-)r], u=UU Q ,0 = 



u, 



0 



u, 



0 



T-T 

T -T 

w oo 



c- c 



, c = 



00 



c — c 



(8) 



Using equations 5,6,7,8 the Esq. 2,3&4 can be written in the following dimensionless form: 



Ps 

H + H-) 



J 



V 



dx dZ 



1 



d 2 U 



{l-tji) 15 dZ 2 



+ 



l-<j> + <j> 



(Pfir)f 



8 cosy + 



l-<j> + <j> 



(pPc) 



(PPc)f 



ccosy - MU 



(9) 



1- <j) + <j) 



{ f_ c _ph 
(pc p ) f 



V 



dr 



dZ 



1 



J 



k nf d 2 9 



k 



f dZ 



1 



Q H 0 



1-0 + 0 



ipc p ) f 



r 8c dc^ 

*3 



v 



dr 



dZ 



J 



Sc 



d 2 c 
dZ 2 



(11) 



(10) 



Where the corresponding boundary conditions (5) can be written in the dimensionless form as: 



U(z,t) = 0, 6{z,t) = 0, c(z,t) =0 for t < 0 V z 



1/(0,0 = U 0 



, £ i int , -ints 
1+— (e +e ) 

2 



, 0(0,0 = 1, c(0,f) = 1 



£/(oo,0->0, 0(oo,O ->0,c(oo,0->0 



Vr>0 



(12) 



Here p r is the Prandtl number, S is the suction (S>0) or injection (S<0) parameter, M is the magnetic 

parameter, Q H (> 0) is the heat source parameter or Q H (< 0) is the heat sink parameter , S c is the Schimdth 
number, which are defined as: 



p r = — — , S= M = J 



a 



■f u 0 Pf Uq - 

Where the velocity characteristic U 0 is defined as 



u 0 = 



80f( T w - T oo)Yf 



nl/3 



The local Nusselt number Nu in dimension less form: 



Qh = 



kf Vo- 



te = 



d 



f 



(13) 



Nu=- 



k 



nf 



0'(P) 



■f 



(14) 



III RESULTS AND DISCUSSIONS 

The governing equations are solved by using Method of lines with the help of Mathematica package. 
The variations of velocity U and temperature 0 are graphically exhibited and the Heat Transfer rate (Nu) is 
exhibited in Table - II for various values of <f> , S, M, OL ,Q H by keeping Pr = 6.2, nt = nil and s = 0.02. The 

effect of various parameters is as follows. 

The increase of solid volume fraction reduces the velocity Fig.l and enhances the temperature Fig.6. 
The thickness of momentum and the thermal boundary layers decreases with increase in </> . From Fig.2&7 the 

momentum and thermal boundary layers decreases for injection or suction. From Fig. 3 & 8 the momentum and 
thermal boundary layers decreasing for heat sink or source Q H . The variations of velocity and temperature with 
magnetic parameter M are depicted in Figs. 4 & 9. The effects of inclination angle a on velocity and temperature 
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are exhibited in Figs. 5 & 10. The increase in inclination reduces the velocity and enhances the temperature. 
From Fig. 11 it is observe that the increase in diffusivity or the decrease in viscosity increases the velocity. The 
same is observed in diffusion with variation of Sc Fig. 12. 




cp = 0 , 0.05 , 
0.15 



1 — I 1 1 — I 1 1 — I — 

0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 



Fig.l Variation of U with 0 




Fig. 3 Variation of U with Q H 




a = tt/6 , n/4, tt/3 



1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 




S=0, 1,1.5 



H 1 h 



0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 



Fig. 2 Variation of U with S 
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Fig. 4 Variation of U with M 
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H 1 h 
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Fig. 5 Variation of U with a 



Fig. 6 Variation of 9 with (j> 
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S = 0,1 ,1.5 



H 1 — I h 



0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 




Fig. 7 Variation of 0 with S 



Fig. 8 Variation of 0 with Q H 




M = 0,5,10 



0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 




a = tt/6 , tt/4 , tt/3 



0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 



Fig. 9 Variation of 0 with M 



Fig. 10 Variation of 0 with a 




-i 1 1 1 — 
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0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 



Fig. 1 1 Variation of U with Sc 



Fig. 12 Variation of C with Sc 
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Nusselt Number: 

Table II 

S = 1; M = 5; cc = \ Sc = 0.6 



* 




G* = io 


0 


4. 5034 


4. 8933 


0.05 


4. 9525 


5. 4462 


0.15 


5. 8815 


6. 6397 


0.2 


6. 3913 


7. 3178 
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